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Abstract. Measurements of a s , the coupling strength of the Strong Interaction between quarks and gluons, 
are summarised and an updated value of the world average of a a (M z o) is derived. Building up on previous 
reviews, special emphasis is laid on the most recent determinations of a s . These are obtained from r-decays, 
from global fits of electroweak precision data and from measurements of the proton structure function F2, 
which are based on perturbative QCD calculations up to O(a^); from hadronic event shapes and jet 
production in e + e~ annihilation, based on O(ots) QCD; from jet production in deep inelastic scattering 
and from X decays, based on 0(« 2 ) QCD; and from heavy quarkonia based on unquenched QCD lattice 
calculations. Applying pragmatic methods to deal with possibly underestimated errors and/or unknown 
correlations, the world average value of a s (M z o) results in 

q s (M z o) = 0.1184 ±0.0007 . 

The measured values of a B (Q 2 ), covering energy scales from Q = M T = 1.78 GeV to 209 GeV, exactly 
follow the energy dependence predicted by QCD and therefore significantly test the concept af Asymptotic 
Freedom. 

PACS. 12.38.Qk experimental tests of QCD 
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1 Introduction 

Quantum Chromodynamics (QCD) is the gauge field the- 
ory of the Strong Interaction pQ. QCD describes the in- 
teraction of quarks through the exchange of massless vec- 
tor gauge bosons, the gluons, using concepts as known 
from Quantum Electrodynamics, QED. QCD, however, 
is more complex than QED because gluons themselves, 
other than photons in QED, carry the quantum-" charge" 
of the Strong Interaction, such that gluons interact with 
each other. 

As a consequence of the gluon self-coupling, QCD im- 
plies that the coupling strength a s , the analogue to the 
fine structure constant a in QED, becomes large at large 
distances or - equivalently - at low momentum transfers p. 
Therefore QCD provides a qualitative reason for the ob- 
servation that quarks do not appear as free particles, but 
only exist as bound states of quarks, forming hadrons like 
protons, neutrons and pions. Hadrons appear to be neu- 
tral w.r.t. the strong quantum charge. The quark statistics 
of all known hadrons, their production cross sections and 
decay times imply that there are three different states of 
the strong charge. 

In loose analogy to the behaviour of optical colours, 
the strong charge and the corresponding new quantum 
number is called "colour charge" . Quarks carry one out of 



three different colour charges, while hadrons are colour- 
less bound states of 3 quarks or 3 antiquarks ("baryons"), 
or of a quark and an anti-quark ("mesons"). Gluons, in 
contrast to photons which do not carry (electrical) charge 
by themselves, have two colour charges. 

QCD does not predict the actual value of a s , however 
it definitely predicts the functional form of the energy 
dependence of a s . While an increasingly large coupling at 
small energy scales leads to the "confinement" of quarks 
and gluons inside hadrons, the coupling becomes small at 
high-energy or short-distance reactions; quarks and glu- 
ons are said to be "asymptotically free" , i.e. a s — > for 
momentum transfers Q — > 00. 

The value of a s , at a given energy or momentum trans- 
fer scakQ Q, must be obtained from experiment. Deter- 
mining a s at a specific energy scale Q is therefore a funda- 
mental measurement, to be compared with measurements 
of the electromagnetic coupling a, of the elementary elec- 
tric charge, or of the gravitational constant. Testing QCD, 
however, requires the measurement of a s over ranges of 
energy scales: one measurement fixes the free parameter, 
while the others test the specific QCD prediction of con- 
finement and of asymptotic freedom. 



1 "Large" distances As correspond to As > 1 fm, 
momentum transfers to Q < 1 GeV/c. 



"low" 



2 Here and in the following, the speed of light and Planck's 
constant are set to unity, c = % = 1, such that energies, mo- 
menta and masses are given in units of GeV. 



2 



Siegfried Bethke: The 2009 World Average of a B 



In this review the current status of measurements of 
a s is summarised. Theoretical basics of QCD and of the 
predicted energy dependence of a s are given in Section 2. 
Actual measurements of a s are presented in Section 3. 
A global summary of these results and a determination 
of the world average value of a s (M z o) are presented in 
Section 4. Section 5 concludes and gives an outlook to 
future requirements and developments. 



2 Theoretical basics 

The concepts of Quantum Chromodynamics are presented 
in a variety of text books and articles, as e.g. [21[31IH|S1 
?J7], such that in the following, only a brief summary of 
the basics of perturbative QCD and the running coupling 
parameter a s will be given. 



2.1 Energy dependence of a s 

With the value of a s known at a specific energy scale 
Q 2 , its energy dependence is given by the renormalisation 
group equation 



Q 



2 da s (Q 2 ) 



dQ 2 



= 0(a s (Q 2 )) 



(1) 



The perturbative expansion of the f3 function is calculated 
to complete 4-loop approximation [S]: 

0(a s (Q 2 )) = - o a 2 (Q 2 ) - /?i« s 3 (Q 2 ) 

- f3 2 a 4 s (Q 2 ) - fcaUQ 2 ) + 0(a 6 s ) , (2) 
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and Nf is the number of active quark flavours at the en- 
ergy scale Q. The numerical constants in equation [3] are 
functions of the group constants Ca = N and Cf — 
(N 2 — 1)/2N, for theories exhibiting SU(N) symmetry. 
For QCD and 577(3), C A = 3 and C F = 4/3. 

A solution of equation [T] in 1-loop approximation, i.e. 
neglecting /?i and higher order terms, is 



« S (Q 2 



a s (n 2 



l + a s (p 2 )/3 \n- 



(4) 



Apart from giving a relation between the values of a s at 
two different energy scales, /x 2 at which a s is assumed to 
be known, and Q 2 being another scale for which a s is be- 
ing predicted, equation |4] also demonstrates the property 



of asymptotic freedom: if Q 2 becomes large and 0o is pos- 
itive, i.e. if Nf < 17, a s (Q 2 ) will asymptotically decrease 
to zero for Q 2 — > oo. 

Likewise, equation [4] indicates that a s (Q 2 ) grows to 
large values and diverges to infinity at small Q 2 : for in- 
stance, with a s (fi 2 = M| ) = 0.12 and for typical val- 
ues of Nf = 2 ... 5, a s (Q 2 ) exceeds unity for Q < 
0(100 MeV. . . 1 GeV). Clearly, this is the region where 
perturbative expansions in a s are not meaningful any- 
more. Therefore energy scales below the order of 1 GeV are 
regarded as the nonperturbative region where confinement 
sets in, and where equations [T] and |4] cannot be applied. 

Including 0i and higher order terms, similar but more 
complicated relations for a s (Q 2 ), as a function of a s (/i 2 ) 

and of In ^ as in equation [4] emerge. They can be solved 
numerically, such that for a given value of a s (/x 2 ), choosing 
a suitable reference scale like the mass of the Z° boson, 
/i = M z o, a s (Q 2 ) can be accurately determined at any 
energy scale Q 2 > 1 GeV 2 . 
With 



yl 2 = 



; l/(/3 asO 2 )) 



a dimensional parameter A is introduced such that equa- 
tion [4] transforms into 



a s (Q 2 ) 



1 



0oHQ 2 /A 2 



(5) 



Hence, the A parameter is technically identical to the en- 
ergy scale Q where a s (Q 2 ) diverges to infinity. To give a 
numerical example, A ea 0.1 GeV for a s (Mz« = 91.2 GeV) 
= 0.12 and N f = 5. 

In complete 4-loop approximation and using the A- 
parametrisation, the running coupling is given [9] by 



a s (Q 2 ) = 
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(6) 



where L = In ^Q 2 /Aj^j. The first line of equation [6] in- 
cludes the 1- and the 2-loop coefficients, the second line is 
the 3-loop and the third and the fourth lines denote the 
4-loop correction, respectively. 

The functional form of a s (Q), in 4-loop approximation 
and for 4 different values of Ajjg, is diplayed in figure [1] 
The slope and dependence on the actual value of Ajj^ 
is especially pronounced at small Q 2 , while at large Q 2 
both the energy dependence and the dependence on Ajjg 
becomes increasingly feeble. 

The relative size of higher order loop corrections and 
the degree of convergence of the perturbative expansion 
of a s is demonstrated in figure [2l where the fractional 
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Fig. 1. The running of a s (Q), according to equation^ in 4- Fig- 2. Fractional difference between the 4-loop and the 1-, 2- 



loop approximation, for different values of A-rrg. 



and 3-loop presentations of a s (Q), for Nf = 5 and Ajj-g chosen 
such that, in each order, a a (M z o) = 0.1184. 



difference in the energy dependence of a s , (ai 4 loop ^ — 
a (n-loop)y a (4,-ioop) ^ £ Qr n= i 2 and 3, is presented. The 
values oi Ajjg were chosen such that a s (M z o) = 0.1184 in 
each order, i.e., Aj^ = 90 MeV (1-loop), A-^ = 231 MeV 
(2-loop), and Aj^ = 213 MeV (3- and 4-loop). Only the 
1-loop approximation shows sizeable differences of up to 
several per cent, in the energy and parameter range cho- 
sen, while the 2- and 3-loop approximations already re- 
produce the energy dependence of the 4-loop prediction 
quite accurately. 

The parametrisation of the running coupling a s (Q 2 ) 
with A instead of a s (/i 2 ) has become a common standard, 
see e.g. [6]. It will also be adopted here. 



2.2 Quark threshold matching 

Physical observables 1Z, when expressed as a function of 
a s , must be continuous when crossing a quark threshold 
where Nf changes by one unit. This determines that A 
actually depends on the number of active quark flavours. 



A will therefore be labelled to indicate these pecu- 

MS 

liarities. Also the slope of the energy dependence and, in 
approximations higher than 2-loop, the value of a s change 
at the quark thresholds: 

Construction of theoretical predictions which consis- 
tently match at a quark flavour threshold leads to match- 
ing conditions for the values of a s above and below that 
threshold [10 . In leading and in next-to-leading order, the 

matching condition is ai Nf — a^ f . In higher orders, 
however, nontrivial matching conditions apply [9l ll0lHT] . 
Formally these are, if the energy evolution of a s is per- 
formed in n th order (or n loops), of order (n — 1). 

The matching scale (j,( N f> can be chosen in terms of the 
(running) mass m q (^i), or of the constant, so-called pole 
mass M q . For both cases, the relevant matching conditions 
are given in [5] . These expressions have a particluarly sim- 



ple form for the choic^fl 



m q {m q ) or ^ N f> 



M„ 



In this review, the latter choice will be used to perform 



3 The results of reference [9] are also valid for other relations 
between and m q or M q , as e.g. m (JV/) = 2M 9 . For 3- 
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3-loop matching at the heavy quark pole masses, in which 
case the matching condition reads, with a — ct& /"K and 



(Nf-1) 



/tt: 



- = 1 + C 2 a 2 



-0.291667 and C 3 



C 3 a* , 
= -5.32389 



(7) 



1) 



where C2 : 
0.26247 [9]. 

The fractional difference of the 4-loop prediction for 
the running a s , using equation[6]with yl^^. - 5 = 213 MeV 
and 3-loop matching at the charm- and bottom-quark 

pole masses, ^c^ -4 ^ = M c = 1.5 GeV and /x^ = 
Mb = 4.7 GeV, and the 4-loop prediction without apply- 
ing matching and with Nf — 5 throughout are illustrated 
in figure O Small discontinuities at the quark thresholds 

can be seen, such that as ^ < as by about 2 per 
mille at the bottom- and about 1 per cent at the charm- 
quark threshold. The corresponding values of Ajj^ are 

yl^ =4) = 296 MeV and A^L =3) = 338 MeV. In addition 

MS .... M S 

to the discontinuities, the matched calculation shows a 
steeper rise towards smaller energies because of the larger 



values of A^—L ^ and A±— L °' . Note that the step func- 

MS M S 

tion of a s is not an effect which can be measured; the 
steps are artifacts of the truncated perturbation theory 
and the requirement that predictions for observables at 
energy scales around the matching point must be consis- 
tent and independent of the two possible choices of (neigh- 
bouring) values of Nf. 



,(JV/=3) 



2.3 Perturbative predictions of physical quantities 

In perturbative QCD, physical quantities 1Z are usually 
given by a power series in a s (/i 2 ), like 

n 

= Pi (R + i?ia s (/i 2 ) + R 2 (Q 2 /n 2 )a 2 (fi 2 ) + ...0) 

where R n are the nth order coefficients of the perturbation 
series and PiRq denotes the lowest-order value of 1Z. R\ 
is the leading order (LO) coefficient, R2 is called the next- 
to-leading order (NLO), R3 is the next-to-next-to-leading 
order (NNLO) and R 4 the N3LO coefficient. 

QCD calculations in NLO perturbation theory are avail- 
able for many observables 1Z in high energy particle re- 
actions like hadronic event shapes, jet production rates, 
scaling violations of structure functions. Calculations in- 
cluding the complete NNLO are available for some totally 
inclusive quantities, like the total hadronic cross section 
in e + e _ — > hadrons, moments and sum rules of struc- 
ture functions in deep inelastic scattering processes, the 
hadronic decay widths of the Z° boson and of the r lepton. 
More recently, NNLO predictions were provided for exclu- 
sive quantities like hadronic event shape distributions and 



loop matching, differences due to the freedom of this choice 
are negligible. 



oc s [matched] 
oc s [unmatched] 




2.5 3 3.5 
Q [GeV] 



Fig. 3. The fractional difference between 4-loop running of 
ct s (Q) with 3-loop quark threshold matching according to equa- 
tions M and [3 with A^—l =5) = 213 MeV and charm- and 

1-1 * MS 

( Nf —4) 

bottom-quark thresholds at the pole masses, fi c = M c = 



1.5 GeV and ^ 



(N f =5) 



M b = 4.7 GeV (full line), and the 



unmatched 4-loop result (dashed line). 



differential jet production rates in e + e _ annihilation |12j . 
and N3LO predictions for the hadronic width of the Z° 
boson and the r lepton [T3] became available. 



A further approach to calculating higher order cor- 
rections is based on the resummation of logarithms which 
arise from soft and collinear singularities in gluon emission 
[14] . Application of resummation techniques and appro- 
priate matching with fixed-order calculations are further 
detailed e.g. in [7]. 
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2.4 Renormalisation 



tegration limits one thus obtains 



In quantum field theories like QCD and QED, physical 
quantities TZ can be expressed by a perturbation series in 
powers of the coupling parameter a s or a, respectively. 
If these couplings are sufficiently small, i.e. if a s <C 1, 
the series may converge sufficiently quickly such that it 
provides a realistic prediction of TZ even if only a limited 
number of perturbative orders will be known. 

In QCD, examples of such quantities are cross sections, 
decay rates, jet production rates or hadronic event shapes. 
Consider TZ being dimensionless and depending on a s and 
on a single energy scale Q. When calculating 1Z as a per- 
turbation series of a pointlike field theory in a s , ultraviolet 
divergencies occur. These divergencies are removed by the 
"renormalisation" of a small set of physical parameters. 
Fixing these parameters at a given scale and absorbing 
this way the ultraviolet divergencies, introduces a second 
but artificial momentum or energy scale /x. As a conse- 
quence of this procedure, 1Z and a s become functions of 
the renormalisation scale fj,. Since 7Z is dimensionless, we 
assume that it only depends on the ratio Q 2 //i 2 and on 
the renormalized coupling a s (fj, 2 ): 

7Z = 7Z(Q 2 /fi 2 ,a s ); a s = a s (n 2 ). 

Because the choice of /i is arbitrary, however, the ac- 
tual value of the experimental observable TZ cannot de- 
pend on [i, so that 



(9) 

where the derivative is multiplied with /x 2 in order to keep 
the expression dimensionless. Equation [5] implies that any 
explicite dependence of TZ on /i must be cancelled by an 
appropriate /x-dependence of a s to all orders. It would 
therefore be natural to identify the renormalisation scale 
with the physical energy scale of the process, p? = Q 2 , 
eliminating the uncomfortable presence of a second and 
unspecified scale. In this case, a s transforms to the "run- 
ning coupling constant" a s (Q 2 ), and the energy depen- 
dence of TZ enters only through the energy dependence of 
a s (Q 2 ). 

The principal independence of a physical observable 
TZ from the choice of the renormalisation scale fi was 
expressed in equation [9] Replacing a s by a s {^ 2 ), using 
equation [TJ and inserting the perturbative expansion of 
TZ (equation [8]) into equation [9] results, for processes with 
constant Pi, in 



= fi 



,ago 



2 " " ' ' i / 2 \ 2 9R\ 2,2, 

1 Mm jm ...... 



dfi 2 



2dR2 R B 



2(>11 '' -[R 1 p 1 +2R 2 p ] 



dfi 2 



(10) 



Rq 


= const. , 
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Solving this relation requires that the coefficients of a™(/i 2 ) 
vanish for each order n. With an appropriate choice of in- 



(11) 



as a solution of equation [TO] 

Invariance of the complete perturbation series against 
the choice of the renormalisation scale /i 2 therefore implies 
that the coefficients R n , except Rq and R\, explicitly de- 
pend on ji 2 . In infinite order, the renormalisation scale 
dependence of a s and of the coefficients R n cancel; in any 
finite (truncated) order, however, the cancellation is not 
perfect, such that all realistic perturbative QCD predic- 
tions include an explicit dependence on the choice of the 
renormalisation scale. 

The scale dependence is most pronounced in leading 
order QCD because R\ does not explicitly depend on /z 
and thus, there is no cancellation of the (logarithmic) scale 
dependence of a s (n 2 ) at all. Only in next-to- leading and 
higher orders, the scale dependence of the coefficients R n , 
for n>2, partly cancels that of a s (fi 2 ). In general, the de- 
gree of cancellation improves with the inclusion of higher 
orders in the perturbation series of TZ. 

Renormalisation scale dependence is often used to test 
and specify uncertainties of theoretical calculations for 
physical observables. In most studies, the central value of 
a s (/i 2 ) is determined or taken for \i equalling the typical 
energy of the underlying scattering reaction, like e.g. p 2 — 
E 2 m in e + e~ annihilation. Changes of the result when 
varying this definition of /i within "reasonable ranges" are 
taken as systematic higher order uncertainties. 

There are several proposals of how to optimise or fix 
the renormalisation scale, see e.g. [TMTOllTTliro] Unfortu- 
nately, there is no common agreement of how to optimise 
the choice of scales or how to define the size of the corre- 
sponding uncertainties. This unfortunate situation should 
be kept in mind when comparing and summarising results 
from different analyses. 

In next-to-leading order, variation of the renormalisa- 
tion scale is sufficient to assess and include theoretical un- 
certainties due to the chosen renormalisation scheme and 
the limited (truncated) perturbation series. In NNLO and 
higher, howerer, both the renormalisation scale and the 
renormalisation scheme should be varied for a complete 
assessment. While it has become customary to include 
renormalisation scale variations when applying theoreti- 
cal predictions, changes of the renormalisation scheme are 
rarely explored. Instead, the so-called "modified minimal 
subtraction scheme" (MS) J9\ is commonly used in most 
analyses, which is also the standard choice in this review. 
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2.5 Nonperturbative methods 

At large distances or low momentum transfers, a s becomes 
large and application of perturbation theory becomes in- 
appropriate. Nonperturbative methods have therefore been 
developed to describe strong interaction processes at low 
energy scales of typically Q 2 < 1 GeV 2 , such as the frag- 
mentation of quarks and gluons into hadrons ("hadroni- 
sation") and the masses and mass splittings of mesons. 

Hadronisation models are used in Monte Carlo ap- 
proaches to describe the transition of quarks and gluons 
into hadrons. They are based on QCD-inspired mecha- 
nisms like the "string fragmentation" 20,21 or "cluster 
fragmentation" [22] , and are usually implemeted, together 
with perturbative QCD shower and/or (N)LO QCD gener- 
ators, in models describing complete hadronic final states 
in high energy particle collisions. Those models contain 
a number of free parameters which must be adjusted in 
order to reproduce experimental data well. They are in- 
dispensable tools not only for detailed QCD studies of 
high energy collision reactions, but are also important to 
assess the resolution and acceptance of large particle de- 
tector systems. 

Power corrections are an analytic approach to approx- 
imate nonperturbative hadronisation effects by means of 
perturbative methods, introducing a universal, non- per- 
turbative parameter 



1 f^ 1 

ato(fii) = — / dk a s (k) 
M/ Jo 



to parametrise the unknown behaviour of a s (Q) below 
a certain infrared matching scale \ii [23] . Power correc- 
tions are regarded as an alternative approach to describe 
hadronisation effects on event shape distributions, instead 
of using phenomenological hadronisation models. 

Lattice Gauge Theory is one of the most developed 
nonperturbative methods (see e.g. [24]) and is used to cal- 
culate, for instance, hadron masses, mass splittings and 
QCD matrix elements. In Lattice QCD, field operators 
are applied on a discrete, 4-dimensional Euclidean space- 
time of hypercubes with side length a. Finite size lattice 
and spacing effects are studied by using increasing lattice 
sizes and decreasing lattice spacing a, hoping to eventu- 
ally approach the continuum limit. With ever increasing 
computing power and refined Monte Carlo methods, these 
calculations significantly matured over time and recently 
provided predictions of the proton (and other hadron) 
masses to better than 2% [25 j5 and determinations of a s 
from quarkonia mass splittings with a precision of better 
than 1 % 126]. 



3 Measurements of a s 

Since almost 30 years, determinations of a s continue to 
be at the forefront of experimental studies and tests of 
QCD. Increasing precision of QCD predictions and meth- 
ods, improved understanding and parametrisation of non- 
perturbative effects, increased data quality and statistics 



and the availabilty of data over large ranges of energy and 
from a large variety of processes have led to an ever in- 
creasing precision and depth of these studies. The timely 
development of a s determinations was documented and 
summarised in a number of summary articles, see e.g. [27 , 
I7ll6ll28]. on which this review is intended to build up. Since 
about the year 2000, the precision of a s determinations 
and the multitude of results from various processes and 
ranges of energies provided experimental proof [71128] of 
the concept of asymptotic freedom. 

In this review, special emphasis is laid on an update 
of the review from 2006 [28] , concentrating on the most 
recent results which are mostly based on further improved 
theoretical predictions and/or experimental precision: 

— perturbative QCD predictions in complete N3LO for 
the hadronic widths of the Z° boson and the r lepton 
are now available, improving further the completeness 
of the perturbative series and providing increased con- 
trol of remaining theoretical uncertainties; 

— improved lattice QCD simulations with vacuum po- 
larisation from u, d and s quarks, updating previous 
determinations of a s and quoting overall uncertainties 
of less than 1%; 

— an improved extraction of a s from radiative decays of 
the T(ls); 

— a combined analysis of non-singlet structure functions 
from deep inelastic scattering data, based on QCD pre- 
dictions complete to N3LO; 

— a combined analysis of inclusive jet cross section mea- 
surements in neutral current deep inelastic scattering 
at high Q 2 ; 

— determinations of a s from hadronic event shapes and 
jet rates in e + e~ annihilation final states, an impor- 
tant and (experimentally) very precise environment, 
based on the new and long awaited QCD predictions 
in complete NNLO QCD. 

These recent results are superior to and thus replace a 
large number of a s determinations published before 2006 
and summarised in [28] . While those previous measure- 
ments still remain valid, only the new results listed above 
are presented and summarised below. 



3.1 a s from r-lepton decays 

Determination of a s from r lepton decays is one of the 
most actively studied fields to measure this basic quan- 
tity. The small effective energy scale, Q — M T = 1.78 GeV, 
small nonperturbative contributions to experimental mea- 
surements of a total inclusive observable, the normalised 
hadronic branching fraction of r lepton decays, 



Rt = 



r(r — ► hadrons v T ) 



rr- 



e v e v T 



(12) 



invariant mass distributions (spectral functions) of hadronic 
final states of r-decays, and the "shrinking error" effect of 
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Fig. 4. Determinations of a s from hadronic r lepton decays 
13 34 ,35 36 ,37 ,38 . The results are all based on the same 
experimental data and on perturbative QCD predictions to 
C(q*), however vary in preference or range of the perturbative 
expansion and inclusion and treatment of nonperturbative cor- 
rections. The vertical line and shaded band show the average 
value and uncertainty used as overall result from r decays in 
this review. 

the QCD energy evolution of a s towards higher energies @ 
provide the means for one of the most precise determina- 
tions of a s {Mio). Theoretically, R T is predicted to be [29] 

R T = N c Sew \V u d\ 2 (1 + S' EW + <5 per t + ^nonpert) ■ (13) 

Here, S E w = 1.0189(6) [30] and S' EW = 0.001(1) [31] are 
electroweak corrections, \V u d\ 2 — 0.97418(27) [6], 6 pelt and 
<5nonport are perturbative and nonperturbative QCD cor- 
rections. Most recently, S pcT t was calculated to complete 
N3LO perturbative order, 0(a^) [13] : it is of similar struc- 
ture as the one for the hadronic branching fraction Rz of 
the Z° boson. Based on the operator product expansion 
(OPE) [32], the nonperturbative corrections are estimated 
to be small [29], <5 n on P ert ~ -0.007 ± 0.004. A comprehen- 
sive review of the physics of hadronic r decays was given 
in [33] . 

Recently, several authors have revisited the determina- 
tion of a s from r decays 13,34,35,36,37,38]. These stud- 
ies are based on data from LEP [39] and - partly- from 
BABAR [40]. They differ, however, in the detailed treat- 
ment and usage of the perturbative QCD expansion of 
R T . In particular, the usage of either fixed order (FOPT) 
or contour improved perturbative expansion (CIPT), and 
differences in the treatment and inclusion of nonpertur- 
bative corrections, leads to systematic differences in the 
central values of a s (M T ), ranging from 0.316 to 0.344, as 
summarised in figure [4] Results based on FOPT turn out 
to be systematically lower than those using CIPT - a trend 
being known for quite some time, and being actively dis- 
puted in the literature, but not finally being solved. 



According to equations 1 and 2, in leading order, 
Aa B (Q 2 ) /a B (Q 2 ) ~ a B (Q 2 ). Therefore, since a B (Q 2 ) decreases 
by about a factor of 3 when running from Q 2 = M 2 to M§, 
the relative error of a s also decreases by about a factor of 3. 



The results shown in figure IU within their assigned to- 
tal uncertainties, are partly incompatible with each other. 
This is especially true if considering that they are based on 
the same data sets. The main reason for these discrepan- 
cies is the usage of either the FOPT or the CIPT pertuba- 
tive expansions. Note that only the result of Baikov et al. 
[13] averages between these two expansions, and assigns an 
overall error which includes the difference between these 
two. 

In view of these differences and for the sake of including 
the apparent span between different perturbative expan- 
sions in the overall error, the range shown as shaded band 
in figure [4] and the corresponding central value is taken as 
the final result from r-decays, leading to 

a s [M T ) = 0.330 ±0.014 . 

Running this value to the Z° rest mass of 91.2 Ge,V using 
the 4-loop solution of the /3-function (equation 6) with 3- 
loop matching at the heavy quark pole masses M c = 1.5 
GeV and M b = 4.7 GeV, results in 

a s (M z o) = 0.1197 ±0.0016 . 

This value will be included in determining the world av- 
erage of a s (M z o) as described in section 4. 

3.2 a s from heavy quarkonia 

Heavy Quarkonia, i.e. meson states consisting of a heavy 
(charm- or bottom-) quark and antiquark, are a classical 
testing ground for QCD, see e.g. [41]. Masses, mass split- 
tings between various states, and decay rates are observ- 
ables which can be measured quite accurately, and which 
can be predicted by QCD based on both perturbation the- 
ory and on lattice calculations. 

3.2.1 a s from radiative Y decays 

Bound states of a bottom quark and antiquark are poten- 
tially very sensitive to the value of a s because the hadronic 
decay proceeds via three gluons, T — ► ggg —> hadrons. 
The lowest order QCD term (i.e. Pi in equation [8]) for 
the hadronic Y decay width already contains a s to the 
3 rd power. The situation is more complicated, however, 
due to relativistic corrections and to the unknown wave 
function of the Y at the origin. 

The wave function and relativistic corrections largely 
cancel out in ratios of decay widths like 

_ r(Y -> 7 gg) 
K -'~ r(Y^ggg) 

which therefore are the classical observables for precise 
determinations of a s . 

In [42], recent CLEO data [43] are used to determine 
a s from radiative decays of the T(1S). The theoretical 
predictions include QCD up to NLO. They are based on 
recent estimates of colour octet operators and avoid any 
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model dependences. The value obtained from this study 
is 

a s (M z o) = 0.119±°:g8i . 

It is compatible with previous results from similar studies, 
see e.g. [7,28] and references quoted therein. It will be 
included in the calculation of the new world average of 
a s (M z o). 

3.2.2 a s from lattice QCD 

Determinations of a s based on lattice QCD calculations 
have become increasingly inclusive and precise in the past, 
including light quarks (u,d and s) in the vacuum polarisa- 
tion and incorporating finer lattice spacing. 

In a recent study by the HPQCD collaboration [26 , 
the QCD parameters - the bare coupling constant and the 
bare quark masses - are tuned to reproduce the measured 
T' - T meson mass difference. The u, d and s quark masses 
are adjusted to give correct values of various light meson 
masses. With these parameters set, there are no other free 
physical parameters, and the simulation is used to provide 
accurate QCD predictions. Nonperturbative values of sev- 
eral short-distance quantities are computed and compared 
to respective perturbative calculations which are given in 
NNLO perturbation theory. From a fit to 22 short distance 
quantities, the value of 

a s (M z o) = 0.1183 ±0.0008 

is finally obtained. The total error incudes finite lattice 
spacing, finite lattice volume, perturbative and extrapola- 
tion uncertainties. This result will be an important ingre- 
dient of the new world average determined in section 4. 

3.3 a s from deep inlastic scattering 

Measurements of scaling violations in deep inelastic lepton- 
nucleon scattering belong to the earliest methods used to 
determine a s . The first significant determinations of a s , 
being based on perturbative QCD prediction in NLO, date 
back to 1979 [B]. 

Today, a large number of results is available from data 
in the energy (Q 2 ) range of a few to several thousand 
GeV 2 , using electron-, muon- and neutrino-beams on var- 
ious fixed target materials, as well as electron-proton or 
positron-proton colliding beams at HERA. In addition to 
scaling violations of structure functions, a s is also deter- 
mined from moments of structure functions, from QCD 
sum rules and — similar as in e + e~ annihilation — from 
hadronic jet production and event shapes. Improved QCD 
predictions as well as new experimental studies provided 
new results from a combined study of world data on struc- 
ture functions, and from jet production at HERA. 

3.3.1 a s from world data on non-singlet structure functions 

Physical processes in lepton-nucleon and in hadron-hadron 
collisions depend on quark- and gluon-densities in the nu- 
cleon. Assuming factorisation between short-distance, hard 



scattering processes which can be calculated using QCD 
perturbation theory, and low-energy or long-range pro- 
cesses which are not accessible by perturbative methods, 
such cross sections are parametrised by a set of structure 
functions Fi [i= 1,2,3). 

The energy dependence of structure functions is given 
by perturbative QCD. A study [45] of the available world 
data on deep inelastic lepton-proton and lepton-deuteron 
scattering provided a determination of the valence quark 
parton densities and of a s in wide ranges of the Bjorken 
scaling variable x and Q 2 . In the non-singlet case, where 
heavy flavour effects are negligibly small, the analysis is 
extended to 4-loop level, i.e. to QCD in N3LO perturba- 
tive expansion. 

The determination of a s to this level results in 

a s (M z o) = 0.1142 ±0.0023 , 

where the total error includes a theoretical uncertainty of 
±0.0008 which is taken from the difference between the 
N3LO and the NNLO result. This value will be included 
in the determination of the world average of a s (M z o). 

As it appears, fits of a s in determinations of parton 
density functions from purely deep inelastic scattering pro- 
cesses result in somewhat smaller values than those which 
also include hadron collider data, see e.g. [JS] obtaining 
a s (M z o) = 0.1171 ± 0.0037. They also are systematically 
smaller than the world average value of a s (M z o), see be- 
low. 

3.3.2 a s from jet production in deep inelastic scattering 
processes 

Measruements of a s from jet production in deep inelastic 
lepton- nucleon scattering at the HERA collider have been 
and continue to be an active field of research. Inclusive as 
well as differential jet production rates were studied in 
the energy range of Q 2 ~ 10 up to 15000 GeV 2 , based 
on similar jet definitions and algorithms as used in e + e~ 
annihilation. 

In a recent summary and combination [47j of precision 
measurements at HERA, values of a s where determined 
from fits of NLO QCD predictions to data of inclusive jet 
cross sections in neutral current deep inelastic scattering 
at high Q 2 [48149] . The overall combined result, 

a s (M z o) = 0.1198 ±0.0032 , 

has a reduced theoretical uncertainty of ±0.0026 (added 
in quadrature to the experimental error of ±0.0019) com- 
pared to previous combinations, due to carefully selected 
ranges of data in Q 2 ranges where theoretical uncertainties 
are minimal [47] . This combined result will be included in 
the determination of the world average of a s (M z <> ) . 

3.4 a s from hadronic event shapes and jet production 
in e + e~ annihilation 

Observables parametrizing hadronic event shapes and jet 
production rates are the classical inputs for a s studies 
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in e + e~ annihilation. The measurements summarised in 
previous reviews [7,50,51,128] were based on QCD predic- 
tions in NLO, which partly included summation of next- 
to-leading logarithms (NLLA) to all orders. As one of 
the most notable and long awaited theoretical improve- 
ments, complete NNLO predictions became available re- 
cently [12j , which are also matched with leading and next- 
to-leading logarithms resummed to all orders [52] (NNLO 
+ NLLA). 

The advancement in theoretical descriptions was in- 
stantly used to determine a s from data of previous e + c~ 
annihilation experiments, from the PETRA and the LEP 
colliders which operated from 1979 to 1986 and from 1989 
to 2000, respectively. The usage of data of past exper- 
iments demonstrates the need to preserve data as well 
as reconstruction-, simulation- and analysis-software for a 
time-span significantly exceeding the usual ~ 5-year pe- 
riod of post-data taking analysis. 

A re-analysis 53 of the ALEPH data from LEP, in 
the cm. energy range from 90 to 206 GeV, based on six 
event shape and jet production observables, results in 

a s (M z o) = 0.1224 ±0.0039 . 

The total error contains an experimental uncertainty of 
0.0013 and a theoretical uncertainty, mainly from hadro- 
nisation and from renormalisation scale dependences, of 
0.0037. This result is obtained using NNLO+NLLA QCD 
predictions; in NNLO alone, the central value is slightly 
higher (0.1240) and the total error is sligthly smaller (0.0032) 
NNLA terms, although they should provide a more com- 
plete perturbation series, tend to introduce somewhat larger 
scale uncertainties [53]. 

Similar results are available from a re-analysis of data 
from the JADE experiment at PETRA [53] , from six event 
shape and jet observables at six cm. energies in the cm. 
energy range from 14 to 46 GeV: 



analysis is a s (Af z o) = 0.1172 ± 0.0021, whereby the error 
includes a theoretical uncertainty of ±0.0017. Although 
this is formally one of the smallest errors quoted on mea- 
surements of a s (M z o), this result is not explicitly included 
in the world average calculated below. The reason for this 
decision is two-fold: First, the LEP thrust data are al- 
ready included in the re-analysis of ALEPH data described 
above. Second, this analysis based on effective field theory, 
although being a highly interesting, alternative approach 
to studies based on standard perturbation theory, is not 
yet in a state of comparable reliability because it is based 
on one event shape observable only, and therefore misses 
an important verification of potential systematic uncer- 
tainties. 



a s (M z o) = 0.1172 ±0.0051 . 

The total error contains an experimental uncertainty of 
0.0020 and a theoretical uncertainty of 0.0046. Also this 
result is obtained using QCD predictions in NNLO±NLLA; 
the value for NNLO alone is a s (M z o) = 0.1212 ± 0.0060. 

Both the NNLO+NLLA results from ALEPH and from 
JADE data are retained for the determination of the new 
world average value of a s (Af z o) in this review. Because 
they are based on data at different cm. energy ranges and 
from two independent experiments, they add valuable and 
independent information not only on the world average, 
but also on the experimental verification of the running of 
a s . These two results of a s (Af z o) are included in figure [5J 
the respective values of a s (Q) , obtained at different values 
of cm. energies, are displayed in figure [5] 

Recently, the event shape observable thrust [55] was 
also studied from LEP data using methods of effective 
field theory [56! . Starting from a factorisation theorem 
in soft-collinear effective theory, the thrust distribution 
is calculated including rcsummation of the next-to-next- 
to-next-to- leading logarithms (N3LL). The result of this 



3.5 a s from electroweak precision data 

The determination of a s from totally inclusive observ- 
ables, like the hadronic width of the r lepton discussed 
above, or the total hadronic decay width of the Z° bo- 
son, are of utmost importance because they lack many 
sources of systematic uncertainties, experimental as well 
as theoretical, which differential distributions like event 
shapes or jet rates suffer from. In this sense, the ratio of 
the hadronic to the leptonic partial decay width, Rz = 
r(Z° -> hadrons)/r{Z Q -> c+e"), is a "gold plated" ob- 
servable. 

Since 1994, the QCD correction to Rz is known in 
NNLO QCD [57]. The measured value from LEP, R z = 
20.767± 0.025 [58], results in a s (M z a) = 0.1226 ± 0.0038, 
where the error is experimental. An additional theoretical 
uncertainty was estimated [7] as 1° 0005 ■ 

As already mentioned, the full N3LO prediction of Rz, 
i.e. in 0(af) perturbative expansion, is now available [13] . 
The negative O(a^) term results in an increase of a s (M z a) 
by 0.0005, such that the actual result from the measured 
value of R z is: a s (M z o) = 0.1231±0.0038 . Defining the re- 
maining theoretical uncertainty as the difference between 
the NNLO and the N3LO result, the theory error would 
not visibly contribute any more, given the current size of 
the experimental error on a s (M z o) of ±0.0038. 

A more precise value of a s (M z o) can be obtained from 
general fits to all existing electro-weak precision data, us- 
ing data from the LEP and the SLC e + c~ colliders as well 
as measurements of the top-quark mass and limits on the 
Higgs boson mass from Tevatron and LEP. Such global 
fits result in values of a s (M z o ) with reduced experimental 
errors. These values, however, were consistently smaller 
than (but still compatible with) the ones obtained from 
Rz alone, see e.g. [58] . 

A recent revision of the global fit to electroweak pre- 
cision data [59], based on a new generic fitting package 
G fitter [60], on the up-to-date QCD corrections in N3LO, 
on proper inclusion of the current limits from direct Higgs- 
searches at LEP and at the Tevatron and on other im- 
proved details, results in 



a s (M z o) = 0.1193; 



-0.0028 
-0.0027 



±0.0005 
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Table 1. Summary of recent measurements of a s (M z o). All eight measurements will be included in determining the world 
average value of a s (M z o ) . The rightmost two columns give the exclusive mean value of a s (M z o ) calculated without that particular 
measurement, and the number of standard deviations between this measurement and the respective exclusive mean, treating 
errors as described in the text. The inclusive average from all listed measurements gives a s (M z o) = 0.11842 ± 0.00067. 



Process 


Q [GeV] 


a B (Q) 


a s (A/ z o) 


excl. mean a s (M z o) 


std. dev. 


r-decays 




1.78 


0.330 ±0.014 


0.1197 ±0.0016 


0.11818 ± 0.00070 


0.9 


DIS [F 2 ] 




2 - 15 




0.1142 ±0.0023 


0.11876 ±0.00123 


1.7 


DIS [e-p -> 


jets] 


6 - 100 




0.1198 ±0.0032 


0.11836 ±0.00069 


0.4 


QQ states 




7.5 


0.1923 ±0.0024 


0.1183 ±0.0008 


0.11862 ± 0.00114 


0.2 


T decays 




9.46 


n 18/I+ 0015 

U - 1 ° 4 -0.014 


n llq +o.oo6 
u.ny_ 005 


0.11841 ± 0.00070 


0.1 


e + e~ [jets k 


z shps] 


14 - 44 




0.1172 ±0.0051 


0.11844 ±0.00076 


0.2 


e + e~ [ew] 




91.2 


0.1193 ±0.0028 


0.1193 ±0.0028 


0.11837 ± 0.00076 


0.3 


e + e~ [jets k 


z shps] 


91 - 208 




0.1224 ±0.0039 


0.11831 ±0.00091 


1.0 



where the first error is experimental and the second is 
theoretical, estimated by the difference of the results in 
NNLO and in N3LO QCD. 



4 The 2009 world average of a s (M z o) 

The new results discussed in the previous section are sum- 
marised in table 1 and in figure Since all of them are 
based on improved theoretical predictions and methods, 
and/or on improved data quality and statistics, they su- 
persede and replace their respective precursor results which 
were summarised in a previous review |28j . While those 
previous results continue to be valid measurements, they 
are not discussed in this review again, and they will not be 
included in the determination of a combined world average 
values of a s (M z o ) , according to the following reasons: 

1. from each class of measurements, only the most ad- 
vanced and complete analyses shall be included in the 
new world average; 

2. older measurements not being complemented or super- 
seded by the most recent results listed above, as e.g. 
results from sum rules, from singlet structure functions 
of deep inelastic scattering, and from jet production 
and b-quark production at hadron colliders, are not 
included because their relatively large overall uncer- 
tainties, in general, will not give them a sizable weight 
but will complicate the definition of the overall error 
of the combined value of a s (M z o); 

3. restricting the new world average to the most recent 
and most complete (i.e. precise) results allows to ex- 
amine the consistency between the newest and the pre- 
vious generations of measurements and reviews. 



4.1 Numerical procedure 

The average x of a set of n different, uncorrelated mea- 
surements xi of a particular quantity x with individual 
errors or uncertainties, Axi is commonly defined using 
the method of least squares (see e.g. [6]: For Xi being 



independent and statistically distributed measures with a 
common expectation value x but with different variances 
Axi, the weighted average is defined by 



EI 



WiX % 



En 
i=l w i 

and the variance Ax of x is minimised by choosing 
1 1 



Ax 



^2 



En l 
i=l "Ax 2 



i.e. Wi 



Ax? 



(14) 



(15) 



The quality of the average is defined by the \ 2 variable 

{Xi -x) 2 



x 2 = 



E 



(16) 



which is, for uncorrelated data, expected to be equal to 
the number of degrees of freedom, n^: 



X 



1 



The results summarised in table 1, however, are not 
independent of each other. They are, in the most general 
sense, correlated to an unknown degree. While the statis- 
tical errors of the data and the experimental systematic 
uncertainties contained in the errors are independent and 
uncorrelated, the theoretical uncertainties are very likely 
to be (partly) correlated between different results, because 
they are all based on applying the same underlying the- 
ory, i.e. QCD, and similar methods to obtain estimates of 
theoretical uncertainties are being used. 

For some observables, like e.g. the hadronic widths of 
the Z° boson and the r lepton, the theoretical predictions 
and hence, their uncertainties, are known to be correlated 
by almost 100%. For other cases, like the results based on 
lattice QCD and those based on QCD perturbation the- 
ory, it can be assumed that their theoretical uncertainties 
are not correlated at all. In addition to the inherent lack of 
knowledge of theoretical correlations, estimates of theoret- 
ical uncertainties, in general, are performed in widely dif- 
ferent ways, using different methods and different ranges 
of parameters. 
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The presence of correlated errors, if using the equations 
given above, is usually signalled by x 2 < n df- Values of 
X 2 > ndf, in most practical cases, are a sign of possibly 
underestimated errors. In this review, both these cases are 
pragmatically handled in the following way: 

In the presence of correlated errors, described by a 
covariance matrix C, the optimal procedure to determine 
the average x is to minimise the % 2 function 



E ( x i — X ){G 1 )ij( x j 



x) 



which leads to 



and 



Ax 



7^2 



The choice of C„ = Axf and dj = for i ^ j re- 
tains the uncorrelated case given above. In the presence 
of correlations, however, the resulting x 2 will be less than 
ndf = n — 1. In order to allow for an unknown common 
degree of a correlation /, the method proposed in [61] will 
be applied by choosing Cij —fx Axi x Axj and adjusting 
/ such that x 2 = n ~ 1- 

For cases where the uncorrelated error determimation 
results in x 2 > n df , and in the absence of knowledge which 
of the errors Axi are possibly underestimated, all individ- 
ual errors are scaled up by a common factor g such that 
the resulting value of x 2 / n df, using the definition for un- 
correlated errors, will equal unity. 

Note that both for values of / > or g > 1, Ax 
increases, compared to the uncorrelated (/ = and g = 1) 
case. 



4.2 Determination of the world average 

The eight different determinations of a s (M z o ) summarised 
and discussed in the previous section are listed in ta- 
ble 1 and are graphically displayed in figure [5] Apply- 
ing equations [TJJ [15] and [T5] to this set of measurements, 
assuming that the errors are not correlated, results in 
an average value of a s (M z o) = 0.11842 ± 0.00063 with 
X 2 /n df = 5.4/7. 

The fact that x 2 < ndf signals a possible correlation 
between all or subsets of the eight input results. Assuming 
an overall correletion factor / and demanding that x 2 — 
ndf = 7 requires / = 0.23, inflating the overall error from 
0.00063 to 0.00089. 

In fact, there are two pairs of results which are known 
to be largely correlated: 

— the two results from e + e~ event shapes based on the 
data from JADE and from ALEPH use the same theo- 
retical predictions and similar hadronisation models to 



T-decays (N3LO) 
Quarkonia (lattice) 
Y decays (NLO) 
DIS F 2 (N3LO) 
DIS jets (NLO) 

e + e~ jets & shps (NNLO) ^ 
electroweak fits (N3LO) 
e + e _ jets & shapes (NNLO) 




0.11 0.12 0.13 

OCs(Mz) 



Fig. 5. Summary of measurements of a s (M z o). The vertical 
line and shaded band mark the final world average value of 
« s (M z o) = 0.1184 ± 0.0007 determined from these measure- 
ments. 



correct these predictions for the transitions of quarks 
and gluons to hadrons. While the experimental errors 
are uncorrelated, the theoretical uncertainties may be 
assumed to be correlated to 100%. The latter accounts 
for about 2/3 to 3/4 of the total errors. An appropriate 
choice of correlation factor between the two may then 
be / = 0.67. 

— the QCD predictions for the hadronic widths of the 
r-lepton and the Z° boson are essentially identical, so 
the respective results on a s are correlated, too. The 
values and total errors of a s (M z o ) from r decays must 
therefore be correlated to a large extend, too. In this 
case, however, the error of one measurement is al- 
most entirely determined by the experimental error 
(Z°-decays) , while the other, from T-decays, is mostly 
theoretical. A suitable choice of the correlation factor 
between both these results may thus be / = 0.5. 

Inserting these two pairs of correlations into the error 
matrix C, the X 2 l n df of the averaging procedure results 
in 6.8/7, and the overall error on the (unchanged) central 
value of a s (M z o) changes from 0.00063 to 0.00067. There- 
fore the new world average value of a s (Mzo) is defined to 
be 



a s (M z o) = 0.1184 ±0.0007. 



For seven out of the eight measurements of a s (M z o), 
the average value of 0.1184 is within one standard devi- 
ation of their assigned errors. One of the measurements, 
from structure functions [45) , deviates from the mean value 
by more than one standard deviation, see figure [5] 

The mean value of a s (Af z o) is potentially dominated 
by the a s result with the smallest overall assigned un- 
certainty, which is the one based on lattice QCD [26]. In 
order to verify this degree of dominance on the average 
result and its error, and to test the compatibility of each 
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of the measurements with the others, exclusive averages, 
leaving out one of the 8 measurements at a time, are cal- 
culated. These are presented in the 5 th column of table 1, 
together with the corresponding number of standard de- 
viations between the exclusive mean and the respective 
single measurement. 

As can be seen, the values of exclusive means vary only 
between a minimum of 0.11818 and a maximum 0.11876. 
Note that in the case of these exclusive means and ac- 
cording to the "rules" of calculating their overall errors, 
in four out of the eight cases small error scaling factors 
of g = 1.06. ..1.08 had to be applied, while in the other 
cases, overall correlation factors of about 0.1, and in one 
case of 0.7, had to be applied to assure x 2 / n df = 1- Most 
notably, the average value a s (M z o) changes to cv s (M z o) = 
0.1186±0.0011 when omitting the result from lattice QCD. 



5 Summary and Discussion 

In this review, new results and measurements of a s are 
summarised, and the world average value of a s (M z o), as 
previously given in [7ll28ll6]. is updated. Based on eight 
recent measurements, which partly use new and improved 
N3LO, NNLO and lattice QCD predictions, the new av- 
erage value is 

a s (M z o) = 0.1184 ±0.0007 , 
which corresponds to 

A^L = (213 ± 9 ) MeV . 

MS v 7 

This result is consistent with the one obtained in the pre- 
viuos review three years ago [28] , which was a s (A/ z o) = 
0.1189±0.0010. The previous and the actual world average 
have been obtained from a non-overlapping set of single 
results; their agreement therefore demonstrates a large de- 
gree of compatibility between the old and the new, largely 
improved set of measurements. 

The individual mesurements, as listed in table 1 and 
displayed in figure [5l show a very satisfactory agreement 
with each other and with the overall average: only one 
out of eight measurements exceeds a deviation from the 
average by more than one standard deviation, and the 
largest deviation between any two out of the eight results, 
namely the ones from r decays and from structure func- 
tions, amounts to 2 standard deviations El- 
There remains, however, an apparent and long-standing 
systematic difference: results from structure functions pre- 
fer smaller values of a s (M z o) than most of the others, i.e. 
those from e + e~ annihilations, from r decays, but also 
those from jet production in deep inelastic scattering. This 
issue apparently remains to be true, although almost all of 
the new results are based on significantly improved QCD 

5 The number of standard deviations is defined as the 
square-root of the value of x 2 - 

6 assuming their assigned total errors to be fully uncorre- 
cted. 



predictions, up to N3LO for structure functions, r and Z° 
hadronic widths, and NNLO for e + e~ event shapes. 

The reliability of "measurements" of a s based on "ex- 
periments" on the lattice have gradually improved over 
the years, too. Including vaccum polarisation of three light 
quark flavours and extended means to understand and cor- 
rect for finite lattice spacing and volume effects, the overall 
error of these results significally decreased over time, while 
the value of a s (M z o ) gradually approached the world aver- 
age. Lattice results today quote the smallest overall error 
on a s (Af z o); it is, however, ensuring to see and note that 
the world average without lattice results is only marginally 
different, while the small size of the total uncertainty on 
the world average is, naturally, largely influenced by the 
lattice result. 




Fig. 6. Summary of measurements of a s as a function of the 
respective energy scale Q. The curves are QCD predictions for 
the combined world average value of a a (M z o), in 4-loop ap- 
proximation and using 3-loop threshold matching at the heavy 
quark pole masses M c = 1.5 GeV and Mb — 4.7 GeV. Full sym- 
bols are results based on N3LO QCD, open circles are based on 
NNLO, open triangles and squares on NLO QCD. The cross- 
filled square is based on lattice QCD. The filled triangle at 
Q = 20 GeV (from DIS structure functions) is calculated from 
the original result which includes data in the energy range from 
Q =2 to 170 GeV. 



In order to demonstrate the agreement of measure- 
ments with the specific energy dependence of ce s predicted 
by QCD, in figure [6] the recent measurements of a s are 
shown as a function of the energy scale Q. For those results 
which are based on several a s determinations at different 
values of energy scales Q, the individual values of a s (Q) 
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are displayed. For the value from structure functions such 
a breakup is not possible; instead, the corresponding re- 
sult derived for a typical energy scale of Q = 20 GeV is 
displayed. 

The measurements significantly prove the validity of 
the concept of asymptotic freedom; they are in perfect 
agreement with the QCD prediction of the running cou- 
pling. This is further corroborated by figure El where a 
selected sample of the measurements is plotted, now as 
a function of 1/logQ, in order to demonstrate the data 
reproducing the specific logarithmic shape of the running 
as predicted by QCD, signalling that indeed a s (Q) — > 
for Q — > oo. 



O ~i ^ r 

0.40 ,mr i 



£ Q [GeV] 



<X S (Q) 



0.30 



0.20 



0.10 



0.00 




QCD 0(<) A'^s = 213 MeV 
no flavour threshold matching 



1 / log(Q/GeV) 



Fig. 7. Selected measurements of a s , as a function of the in- 
verse logarithm of the energy scale Q, in order to demonstrate 
concordance with Asymptotic Freedom. The full line is the 
QCD prediction in 4-loop approximation with 3-loop threshold 
matching at the heavy quark pole masses. The dashed line indi- 
cates extrapolation of the 5-flavour prediction without thresh- 
old matching. 



What are the future prospects of measurements of a s ? 
With the given degree of data and theory precision, fur- 
ther improvements will be difficult and may take quite 
some time. Experimentally, a linear e + e~ collider, espe- 
cially if run in the "Giga-Z" mode, has the potential to 
decrease the dominating experimental error of a s (Mio) 



from the measurement of Rz, down and below its theo- 
retical uncertainty which currently is assumed to be, in 
N3LO, ±0.0005. 

While it is unlikely that QCD perturbation theory will 
improve to yet one order higher than the existing N3LO or 
NNLO predictions, improvements are likely, and actually 
are very neccessary, for QCD predictions of jet produc- 
tion in deep inelastic scattering and in hadron collisions, 
where calculations currently are limited to NLO. The pre- 
cision of QCD tests, but also the sensitivity for observing 
new physics signals at the LHC, will largely depend on a 
further advancement of QCD predictions for hadron colli- 
sions. 

Future improvement of theoretical predictions and mod- 
els require the conservation of data and of reconstruction 
and simulation code of current and past experiments; es- 
pecially in the case of deep inelastic scattering data, re- 
application of improved predictions and models carry a 
large potential for future advancements in this field. 

Last but not least, further developments of non- per- 
turbative methods are mandatory to bridge the gap be- 
tween quarks and gluons and their final states, hadrons. 
They may in fact shed more light into the systematic 
differences between some classes of measurements as dis- 
cussed above. 
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